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ABSTRACT 

We report the first investigation of cool-core properties of galaxy clusters selected via their Sunyaev- 
Zel'dovich (SZ) effect. We use 13 galaxy clusters uniformly selected from 178 deg 2 observed with the 
South Pole Telescope (SPT) and followed-up by the Chandra X-ray Observatory. They form an 
approximately mass-limited sample (> 3 x 10 14 M Q /i7 ) spanning redshifts 0.3 < z < 1.1. Using 
previously published X-ray-selected cluster samples, we compare two proxies of cool-core strength: 
surface brightness concentration (csb) and cuspiness (a). We find csb is better constrained. We 
measure csb for the SPT sample and find several new z > 0.5 cool-core clusters, including two 
strong cool cores. This rules out the hypothesis that there are no z > 0.5 clusters that qualify as 
strong cool cores at the hAo level. The fraction of strong cool-core clusters in the SPT sample in 
this redshift regime is between 7% and 56% (95% confidence). Although the SPT selection function 
is significantly different from the X-ray samples, the high-z csb distribution for the SPT sample is 
statistically consistent with that of X-ray-selected samples at both low and high redshifts. The cool- 
core strength is inversely correlated with the offset between the brightest cluster galaxy and the X-ray 
centroid, providing evidence that the dynamical state affects the cool-core strength of the cluster. 
Larger SZ-selected samples will be crucial in understanding the evolution of cluster cool cores over 
cosmic time. 

Subject headings: cooling flows - galaxies: clusters: general - X-rays: galaxies: clusters 
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1. INTRODUCTION 

Galaxy clusters grow over cosmic time through mergers 
with other galaxy clusters as well as through the accre- 
tion of gas and individual galaxies from the surrounding 
environment. On timescales of a few Gyr, radiative cool- 
ing due to X-ray emission from the intracluster medium 
(ICM) w ould give rise to a "coo l ing flow" to the clus 



ter co re (Fabian & Nulsen 1977 Mathews & Bregman 
19781, if it were not countered by a heating mechanism. 



These cooling flows are not observed; instead, the cores of 
some clusters are found to undergo only moderate cooling 
( Kaastra et al.||2001 Peterson et al.|[200T Tamura et al. 
2U01D. Such galaxy clusters are known as "cool-core" 



clusters ( Molendi fe Pizzolato||200T ). Other clusters ex- 
hibit little to no cooling m their core (i.e., noncool-core 
clusters). These cooling properties suggest that there 
must be processes in every cluster that are strong enough 
to either regulate cooling flows or completely prevent 
them. Such processes are not fully understood and it 
is still uncertain how they evolve and affect cluster for- 
mation over time. 

The important astrophysical processes that counteract 
cool-core formation are typically thought to fall under 
three broad categories: feedback from active galactic nu- 
clei (AGN), preheating of the cluster gas, and cluster 
mergers. AGN feedback in the cluster's brightest cluster 
galaxy (BCG) has been shown to be capa ble of regulat 



ing cooling flows in cool-co re clusters (see Fabian 2012 
McNamara & Nulsen 2012 for reviews). In some cases 
the feedback may b e strong enough to d isrupt a cool 
core completely (e.g. McDonald et al.|2011 ), though this 
phenomenon is likely limited to lower mass galaxy clus- 
ters and groups. Additionally, AGN feedback may drive 
turbulence in the ICM. This has been shown to suppress 
heat-flux-driven buoyancy instabilities, resulting in effec- 
tive transfer o f heat from the oute r radii and disrupting 
the cool core (Parrish et al. 2012 1. Heating of the intr- 
acluster gas during early stages of the cluster has been 
shown to affect the formation of cool cores as well (e.g. 
Kaiser||199l| |McCarthy et al.||2008[ |Sun||2009| ). Cluster 
mergers can also di srupt cool cores byshock-heating and 



turbu lent mixing (Leccardi et al.| 



2010 



DIO] i Rossetti et al., 
201 1| , a process that has been reproduced in simulations 
e^lMcGlynn fc Fabian|1984l|G6mez et al.|2002||ZuHone| 



et ah] 2TJTD ). Whether a, merger can destroy a cool core 
likely depends on the strength of the cool core, the mass 
ratio of the merging clusters, and the geometry of the 
impact. 

Studying the evolution of clusters can provide insight 
regarding the relative importance of these processes in 
cool-core formation. Given that cool cores develop over 
a central cooling time of typically a few Gyr, one ex- 
pects there to be fewer cool cores at times closer to 
the epoch of galaxy cluster formation. Simulations pre- 
dict a significantly high er cluster merger rate in the past 
(Gottlober et al. 2001). If mergers play a strong role 
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in the disruption of cool-core galaxy clusters, this too 
suggests that the fraction of galaxy clusters with cool 
cores should be lower at high redshifts than in local sam- 
ples. Indeed, studies of the evolution of the cool-core 
fraction find a much lower fraction at z = 0.5 (McDon- 



z > 0.5 JVikhlinin et al. 2007 


Santos et al.(2010 


Samuele 


et al. 


201ip. To date, only a small number o 


[' z > 0.5 



report ed (e.g. Siemiginowska et al. 2010 |Russell et al. 



20121, with the most dramatic, confirmed strong cool 
core coming from the Sout h Pole Telescope (SPT) sur 



vey ( |McDonald et al.|2012[ ). 

Understanding cool-core evolution is complicated by 
selection biases of cluster samples. One can generically 
expect an X-ray-selected sample to be biased toward 



selecting cool-core galaxy clusters (Hudson et al. 2010 
Eckert et al. 2011 ) due to their higher A.-ray surface 
brightness and luminosity as compared to a noncool- 



core galaxy cluster of the same mass (e.g. O'Hara et al. 
2006|). However, there are competing effects due to X 



ray emission from AGNs, whi ch are expected to be more 



prevalent at higher redshifts ( [Russell et al.||2012[ ). The 

;he 



bias may be complicated further by the ways m which 
different cluster-finding algorithms differentiate between 
point sources (e.g., AGNs, X-ray binaries) and extended 
sources (e.g., nearby galaxies, groups, and clusters). Fur- 
thermore, the classification of the cool-core strength of a 
cluster varies between surveys of different angular reso- 
lution and the method used to characterize the cool core. 

Given the complex effects associated with X-ray selec- 
tion, it is important to investigate the cool-core fraction 
and its evolution using an independent selection method. 
In this paper, we study the cool-core properties of galaxy 
clust ers selected from their Suny aev-Zerdovich (SZ) ef- 



fect ( |Sunyaev fc Zel'dovichl[l972| ) signature in the SPT 
survey. At z > 0.3, this SPT selection is nearl y redshift 



independen t and nearly constant in mass (e.g. Reichardt 
et al. 2012). The SZ effect selection is expected t o be 
relatively insensitive to non-gravitational physics (|Na- 



gai|2006|) , the dynamical state of clusters (|Jeltema et al 



20081), radio contamination from point sources and B(JGs 



(|Lin et ah 



et al.||2005 



to 



2009), and the presence of cool cores (Motl 
Therefore, an SZ-cluster survey is expected 
Be a useful tool to study the redshift evolution of 
galaxy cluster properties. This work provides the first 
glimpse of the cool-core properties of a sample of galaxy 
clusters selected from the SZ effect. 

The layout of this paper is as follows. In Section [2] we 
detail the observations used and describe the data reduc- 
tion procedures. In Section [3J we present the steps used 
to make our measurements as well as compare two meth- 
ods used to characterize cool-core strengths that are suit- 
able for high- redshift clusters. In Section |4j we present 
the results of our measurements and investigate the im- 
plications for the cool-core fraction at high redshifts. In 
Section [5j we investigate the relationship between a clus- 
ter's coohcore strength and the offset of its BCG. Finally, 
in Section [6j we conclude our analyses and present future 
studies and applications. 
In this analysis, we assume the be st-fit 

WMAP7+BAO+ Hn flat ACDM cosmology fiio- 

matsu et al. 2011 ) with Hubble parameter 
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H = 70 .4 km s" 1 Mpc" 1 , matter density fl M = 0.272, 
and dark energy density — 0.728. 

2. OBSERVATIONS 

We analyze Chandra data for two different cluster 
sampl es, an SZ-selected sample from |Andersson et ah 



Table 1 

Chandra Observation IDs 



(2011 hereafter, All) an d an X-ray-selected sample 



from | Vikhlinin et al. (2009 hereafter, V09). In this sec- 
tion we describe each cluster sample and the Chandra 
X-ray data reduction. 

2.1. Cluster Samples 

We first describe the SZ-selected sample from All. 
This cluster sa mple is a subset of the SP T cluster cata- 
log described by Vanderlinde et al. ( 2010[ ), which consists 
of 21 clusters selected by their SZ-significance (£) from 
178 deg 2 of sky surveyed by the SPT. From the subsam- 
ple of 17 clusters at z > 0.3 and £ > 5.45, 15 clusters 
were selected for X-ray ob servations, the results of which 
are discussed in All and Benson et al. (2011). Of the 
15 observations, 13 were carried out by Chandra and two 
were carried out by XMM-Newton. For this study, we 
use the 13 clusters with Chandra observations because 
only Chandra provides the spatial resolution needed for 
our analysis. These 13 clusters form a nearly complete 
mass- limited sample (called the "SPT sample"). 

In addition to these, we analyze SPT-CL J2106-5844, 
a massive gal axy cluster at z = 1.13 discovered by the 
SPT survey (Foley et al. 20111. However, this cluster 
is not included in analyses invo lving the distribution of 
cool-core strengths (Section |4.3[ ) as it is not a member of 
the mass-limited dataset. 

We also analyze 41 galaxy clusters, based largely from 
the high-redshift sample of the Chandra Cluster Cos- 
mology Project 32 , known hereafter as the CCCP high-z 
sample (V09). This X-ray-selecte d sample is the sub- 
set of clusters in the 400d survey ( |Burenin et al. 20071 
at z > 0.35 and above a redshitt-dependent X-ray- 
luminosity threshold. We include five additional clusters 
from the 400d survey with Chandra data available not 
presented in V09. 33 



2.2. Chandra Data Reduction 

The Chandra observations used in this study are listed 
m Table [1] Data are reduced using the Chandra soft- 
ware version CIA0 4.4 and CALDB 4.4.8. For all ob- 
servations, Level = 2 event files arc generated with the 
chandra_repro script. Exposure corrections are applied 
using f luximage in the 0.5 — 2 keV band with exposure 
maps calculated at 1.5 keV. We employ a two-step pro- 
cedure in removing point sources. First, candidate point 
sources are identified using wavdetect. The results are 
then visually inspected and false detections within the 
cluster are ignored. Using dmfilth, proper detections 
are replaced with a level determined from an elliptical 
annulus centered on the point source. Background levels 
are determined in each observation from several regions 
located on the same chip as the cluster source without 
point source detections. The regions are positioned far 

32 http: / /hea- www. harvard.edu/400d/cosm/ 

33 These clusters are included to increase the samp le size as well 
as to remain consistent with the dataset presented in Santos et al. 
I |2010| . 



Cluster 


Observation ID 


SF1-CL J0000— 5748 


9335 


SF1-CL J0509— 5342 


9432 


SPT-CL J0516-5430 


9331 


SPT-CL J0528-5300 


12092, 10862, 11996, 9341, 11874 


SPT-CL J0533-5005 


11748, 12001 


SPT-CL J0546-5345 


9336, 9332, 10864, 10851 


SPT-CL J0551-5709 


11871 


SPT-CL J2106-5844 


12180 


SPT-CL J2331-5051 


11738, 9333 


SPT-CL J2337-5942 


11859 


SPT-CL J2341-5119 


11799, 9345 


SPT-CL J2342-5411 


11870, 11741, 12014 


SPT-CL J2355-5056 


11998, 11746 


SPT-CL J2359-5009 


9334, 11742, 11864 


CL J0340— 2823 


5780 


CL JUoUz— U42o 


5782 


CL M2\2-\-276o 


5767 


CL J0350— 3801 


7227 


CL J0318— 0302 


5775 


CL J1514+363b 


800 


CL J0159+0030 


5777 


CL J0958+4702 


5779 


/~\ t inoAn i iioi i 

CL J0809+2811 


5774 


/~1 T TT A 1 a i A A AC 

uL J141o+444o 


541 


CL J1312+3900 


5781 


CL JlUUo+oZOO 


5776 


/~i t Tm a i ono a 

CL JU141 — o(Jo4 


5778 


/~1 T T -i r-7C\ 1 i C A 1 A 

CL J17U1+D414 


547 


/~1 T T -I C A -1 1 (rim 

CL Jlo41+4UUl 


3575 


CL JUb22— ot)24 


492b, 583/ 


/~i t t i oo o i n'7nn 

CL J1222+27U9 


5766 


CL J0355— 3741 


5761 


CL JUo5o+o7oy 


4925, 5765 


CL JUooo— 24oo 


5764 


CL J0926+1242 


4929, 5838 


CL J0030+2bl8 


5762 


CL J10U2+D858 


5773 


CL J 1524+0957 


1664 


CL J1357+6232 


5763, 7267 


CL Jloo4— U221 


4yo2, OooO 


CL J 1117+1744 


4933, 583b 


/~i T t i i on i o o o /"> 

CL J1120+2326 


3235 


^~ "T Tfl01 £i 1 r 7A r 7 

CL J021D— 1/4/ 


t^.*7iiC\ /?QOQ 

5/60, b39o 


CL J0521-2530 


5758 


CL J0956+4107 


5294, 5759 


CL J0328-2140 


5755, b258 


CL J1120+4318 


5771 


CL J1334+5031 


5772 


CL J0542-4100 


914 


CL J1202+5751 


5757 


CL J0405-4100 


7191 


CL J1221+4918 


1662 


CL J0230+1836 


5754 


CL J0152-1358 


913 


CL J1226+3332 


3180, 5014 



enough from the cluster emission so as to contain negli- 
gible cluster photons and are large enough to adequately 
sample the background level. 

The results enable measurements to be made on the 
background-subtracted, exposure-corrected, flux images 
in the 0.5—2 keV band where each pixel is in units of pho- 
tons cm -2 s . Count-rate errors are determined based 
on Poisson statistics and propagated in the standard way. 

3. MEASURING COOL CORES WITH IMAGING DATA 



Hudson et al. 



(2010) compare several X-ray estima- 
tors of cool-core strength that are applied to a cluster 
sample with a range of X-ray data quali ty and redshifts. 
They find that the X-ray cuspiness (a; |Vikhlinin et al. 
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2007|) and the su rface brightness concentration (csb ; San- 
tos et alJ 2008[ ) are promising cool-core estimators for 
shift ch 



high-redshitt clusters with observations containing rel 
atively few X-ray photons. In the low-redshift regime, 
they find the central cooling time, i coo i to be the best 
cool-core estimator based on the strength of its bimodal- 
ity. Both csb and a have be en shown to correlate well 
with the ce ntral cooling time ( Santos et al.|2008 Hudson 
et al. 20101. In this section, we compare the a and csb 



Table 2 

CCCP High-2 Galaxy Cluster Measurements 



parameters to the central cooling time for a wider sam- 
ple of clusters. We also calculate csb for the portion of 
the CCCP high-z sample not previously published (i.e., 
clusters at z < 0.5). 



Pg 



3.1. Calculating a and csb 
The cuspiness is defined as the slope of the gas density 

(1 lOg Pg 



dlogr 



where the function is evaluated at radius r — 0.04r 50 o, 
^500 being the radius at which the mean density of the 
enclosed mass is 500 ti mes that of the critic al density 
at the object's redshift ( |Vikhlinin et al.||2007| . This ra- 
dius is close enough to the cluster core to sample the 
areas of strongest cooling, while still being far enough 
to avoid any flattening of the density profile caused by 
feedback from a central AGN. Calculations of a for the 
CCCP high-z sample are derived from the X-ray surface 
brightness fits and X-ray centers used in V09. 

The csb is defined as the ratio of the soft X-ray flux, 
F, within the inner 40 kpc to the inner 400 kpc 



F, 



csb 



■r<40 kpc 



F, 



r<400 kpc 



( Santos et al.| 2008). These radii were chosen because 
they provide the largest separation of csb values between 
cool-core and noncool-core clusters. Previous studies 
have used either the 0.5 — 5.0 keV energy band or the 
0.5 — 2.0 keV energy ba nd for determining csb- We fol- 
low Santos et al. (2010) in using the 0.5 — 2.0 keV band 



for our csb measurements. 

T o calculate csb in the clusters described in Section 
|2.1| we first estimate the centroid of the X-ray emission 
for each galaxy cluster. This is determined by iterating 
the centroid within an 80 pixel (~40") radius, initially 
centered on the approximate location of the centroid. For 
each iteration, the centroid is calculated with the image 
weighted by r -1 / 2 , where r is radial distance to the cen- 
ter of the previous iteration. For clusters with multiple 
observations, an inverse variance weighted final csb is de- 
termined using the csb values measured in the individ- 
ual observations. In these clusters, the X-ray centroid is 
taken to be the mean centroid of all observations, with- 
out any weighting applied. 

When comparing X-ray measurements across different 
redshift regimes, a i^-correction is commonly applied to 
account for the redshift dependence of the flux in a given 
ban d. This is generally a small effect for the csb parame- 
ter (Santos et al. 2010). We assume it's negligible for this 
work because we are primarily interested in comparing 
several galaxy cluster samples at similarly high redshifts. 



Cluster 


z 


CSB 


a 




(1) 







CL J0340-2823 a 
CL J0302-0423 
CL J1212+2733 
CL J0350-3801 
CL J0318-0302 
CL J1514+3636 a 
CL J0159+0030 
CL J0958+4702 
CL J0809+2811 
CL J1416+4446 
CL J1312+3900 
CL J1003+3253 
CL J0141-3034 
CL J1701+6414 
CL J1641+4001 
CL J0522-3624 
CL J1222+2709 
CL J0355-3741 
CL J0853+5759 
CL J0333-2456 
CL J0926+1242 
CL J0030+2618 
CL J1002+6858 
CL J1524+0957 
CL J1357+6232 
CL J1354-0221 
CL J1117+1744 a 
CL J1120+2326 
CL J0216-1747 a 
CL J0521-2530 a 
CL J0956+4107 
CL J0328-2140 
CL J1120+4318 
CL J1334+5031 
CL J0542-4100 
CL J1202+5751 
CL J0405-4100 
CL J1221+4918 
CL J0230+1836 
CL J0152-1358 
CL J1226+3332 



0.35 
0.35 
0.35 
0.36 
0.37 
0.37 
0.39 
0.39 
0.40 
0.40 
0.40 
0.42 
0.44 
0.45 
0.46 
0.47 
0.47 
0.47 
0.48 
0.48 
0.49 
0.50 
0.50 
0.52 
0.53 
0.55 
0.55 
0.56 
0.58 
0.58 
0.59 
0.59 
0.60 
0.62 
0.64 
0.68 
0.69 
0.70 
0.80 
0.83 
0.89 



0.114 ± 
0.374 ± 
0.036 ± 
0.073 ± 
0.039 ± 
0.276 ± 
0.175 ± 
0.144 ± 
0.028 ± 
0.149 ± 
0.046 ± 
0.216 ± 
0.058 ± 
0.155 ± 
0.087 ± 
0.048 ± 
0.115 ± 
0.096 ± 
0.025 ± 
0.035 ± 
0.092 ± 
0.040 ± 
0.060 ± 
0.032 ± 
0.054 ± 
0.043 ± 
0.041 ± 
0.027 ± 
0.055 ± 
0.046 ± 
0.040 ± 
0.062 ± 
0.063 ± 
0.068 ± 
0.043 ± 
0.030 ± 
0.073 ± 
0.026 ± 
0.036 ± 
0.027 ± 
0.086 ± 



0.012 

0.024 

0.005 

0.012 

0.007 

0.014 

0.017 

0.016 

0.006 

0.012 

0.009 

0.022 

0.014 

0.013 

0.011 

0.008 

0.014 

0.013 

0.007 

0.007 

0.010 

0.011 b 

0.012 b 

0.006 b 

0.010 b 

0.009 b 

0.010 b 

0.011 b 

0.014 b 

0.007 b 

0.007 b 

0.009 b 

0.005 b 

0.017 b 

0.007 b 

0.008 b 

0.009 b 

0.006 b 

0.009 b 

0.008 b 

0.007 b 



0.920 ± 
1.349 ± 
0.389 ± 
0.126 ± 
0.078 ± 
1.185 ± 
1.102 ± 
0.751 ± 
0.021 ± 
1.008 ± 
0.026 ± 
1.369 ± 
0.563 ± 
1.189 ± 
0.695 ± 
0.671 ± 
0.743 ± 
0.877 ± 
0.190 ± 
0.395 ± 
0.744 ± 
0.358 ± 
0.185 ± 
0.056 ± 
0.632 ± 
0.035 ± 
0.322 ± 
0.280 ± 
0.428 ± 
0.500 ± 
0.026 ± 
0.513 ± 
0.679 ± 
0.381 ± 
0.454 ± 
0.009 ± 
0.293 ± 
0.049 ± 
0.161 ± 
0.102 ± 
0.333 ± 



0.068 
0.053 
0.133 
0.176 
0.102 
0.049 
0.071 
0.148 
0.089 
0.064 
0.036 
0.188 
0.284 
0.043 
0.314 
0.102 
0.132 
0.076 
0.208 
0.173 
0.116 
0.155 
0.147 
0.698 
0.113 
0.434 
0.221 
0.127 
0.203 
0.225 
0.007 
0.185 
0.101 
0.228 
0.136 
0.121 
0.111 
0.032 
1.171 
0.364 
0.075 



Note. — (1) Redshifts from Vikhlinin et al?]j2009[ . 
a Cluster not a member of the original CCCP high-£ set pre- 
sented in V09 

b csb values frompantos et al.|<|201o|. 



Using the above analysis, we report a and csb for the 
CCCP high-z sampl e in Table [2} with cs b values for clus- 



ters at z > 0.5 from Santos et al. (2010b. 



In the following analyses, we adopt the three different 
cool-core regimes defined previously for th ese parameters 
( |Vikhlinin et al.||2007l [Santos et al.||2008| ): the noncool- 
core regime (csb < 0.075; a < 0.5), the moderate regime 
(0.075 < c SB < 0.155; 0.5 < a < 0.7), and the strong 
cool-core regime (csb > 0.155; a > 0.7). 

3.2. Comparison of csb oxid a 

We study the performances of a and csb by first re- 
lating the two parameters to t coo \ in a sample of low- 
redshift galaxy clusters. In Figure Jl] we compa re the 



from ISantos e t al 



;2008) 



measurements of csb and t coo i 

with those of a and t co o i from |Hudson et ai JpjTOt TWe 
note that Santos et al. ( 2008 ) calculate csb in a different 
band of 0.5 — 5.0 keV, however this should not qualita- 
tively change our conclusions because the bulk of ICM 
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10 15 20 

*cooi (Gyr) 

Figure 1. Comparison of a and csb with t r r.r.i ■ Galaxy clu sters 
with a and t coo i values are taken from | Hudson et al,| || 2010j| and 
galaxy clusters with csb and t coo i values are taken from |aantos| 
|et al.| l [2008l ). 



0.5, the median S/N in a is 2.9. However, for all clusters 
in this sample, the median S/N in csb is 7.8. 

As shown in Figure [2j the two parameters do not al- 
ways provide the same classification for a cluster's cool- 
core strength. The classification t hresho lds fo r crr and 
a were defined in |Santos et al. ( 2008 ) and Vikhlinin 
et al. (2007), respectively. The definitions were based 



on separate datasets using characteristics of each param- 
eters' distribution. Therefore, one would not necessarily 
expect perfect agreement between the two parameters. 
Regardless, using these classifications, all strong cool- 
core clusters identified by Csb are also classified as such 
when using a. Therefore, the csb parameter provides a 
more conservative threshold for identifying strong cool- 
core clusters. 

Throughout the rest of the paper, we will use the csb 
parameter to characterize the cool-core strength of both 
the X-ray and SZ-selected samples. The X-ray obser- 
vations of these samples have comparable S/N to the 
samples used in this Section. Therefore, we expect sim- 
ilar performances from a and csb as was shown here, 
where the results suggest that csb is a more robust and 
comparably accurate cool-core proxy. 
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Figure 2. Correlation between a and csr for the CCC P high- 
z clusters. Values for a are from | Vikhlinin et a.l.| (|2007| >. Values 
for csb with z < 0.5 ar e reported in t his st udy, values for csb 
with z > 0.5 are from |Santos et al.| ( |2010fr . The dashed lines 
correspond to the boundaries between noncool cores and moderate 
cool cores. The dash-dotted lines correspond to the boundaries 
between moderate cool cores and strong cool cores. 



emission is in the 0.5 — 2.0 keV band. We find a and 
icooi to have a Spearman's rank correlation coefficient of 
p a = —0.88. The associated p-value corresponds to a 
probability of 4.4 x 10~ 22 that there is no correlation be- 
tween the two parameters. The results are similar for csb 
and tcooi; their Spearman's rank correlation coefficient is 
p CSB — —0.84 with a p- value of 6.3 x 10~ 8 . 

We compare the two parameters directly for the galaxy 
clusters in the CCCP high-z set. As shown in Figure [2j 
the a parameter exhibits larger fractional errors than 
csb in the CCCP high-z sample. This is particularly 
prevalent for low a, indicating that the parameter is less 
constrained in the noncool-core regime. Even excluding 
the 6 clusters where the signal-to-noise ratio is less than 



4. EVOLUTION OF THE COOL-CORE FRACTION 

In this Section, we compare the csb distribution as a 
proxy for the cool-core fraction for various comparable 
cluster samples. We first describe the X-ray-selected 
cluster samples previously studied and then go on to 
present the results for the SPT sample. We compare 
the csb distribution of the X-ray-selcctcd sample with 
that of the SPT sample in different redshift regimes. Fi- 
nally, we discuss systematics that could affect the SPT 
sample. 

4.1. X-ray-Selected Cluster Samples 



Previously, Santos et al. (2010) compared the surface 
brightness concentrations for 57 galaxy clusters. They 
compared both low- and high-rcdshift samples as well 
as clusters selected by different X-ray selection methods. 
All of the clusters were first detected in ROSAT PSPC 
observations and were later followed up with Chandra. 
About half of the clusters (26) are from the 0.05 < z < 
0.22 portion of the CCCP low-z set (V09), a flux-limited 
cl uster sample with m any of the same members presented 
in |Edge et al.| (1990). The other half is comprised of 
clusters in the CCCP high-z set (V09) with z > 0.5 (20 
clusters) 34 and the z > 0.6 clusters fro m the RDCS 35 
( Rosati et al||l998| and the WARPS 36 ( [Horner et al] 
2008[ ) samples that have Chandra observations (15 Clus- 
tersT One of these clusters, WARP J 1415. 1+3612, has 
an updated csb in Santos et al. (2012). Due to the low 
number of counts and similar selection properties in the 
RDCS and WARPS surveys, [Santos et ah] ( |2010[) group 
them together in their analysis as RDCS+\VARPS. For 
consistency, we follow the same practice here. Note that 

34 The 20 clusters from the C CCP high-z sample with z > 0.5 
presented in S antos et al.| ( |20 L0| ) include three clusters not in the 
original CCOP iugh-z sample, as their Chandra observations re- 
vealed they did not meet the flux criteria for membership. 

35 ROSAT Deep Cluster Survey 

36 Wide Angle ROSAT Pointed S urvey, 
http: / /asd. gsfc.nasa.gov/Donald.Horner/warps/index. html 
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four clusters are members of both RDCS+WARPS and 
the CCCP high-z samples. There is one additional high- 
redshift cluster that has a published c sb value, XMMU 



J1230.3+1339 dFassbender et al. 



Santos et al.| ( |20lU ) compare 



2011) 



he csb distributions of 
the CCCP high-"FaEc! RDCS+WARPS samples using a 
K-S test. The result yields a probability of only 0.6% 
for the null hypothesis that the two distributions are 



drawn from the same parent distribution. Santos et al. 
(2010) argue this is due to a bias of the detection al 



gorithm used in CCCP against compact clusters with a 
relati vely high mean surf ace brightness. This conflicts 



with Burenin et al. (2007), who tested the CCCP selec- 
tion algorithm on a morphologically diverse set of cluster 
images, reprojected to redshifts between 0.35 and 0.80. 
The results showed comparable selection efficiencies for 
all of the cluster morphologies, indicating the selection 
method does not present a bias with respect to the cool- 
core strength of the cluster. 

4.2. The SPT Sample 



The SPT sample is described in Section 2.1 It con- 
sists of 13 clusters with X-ray observat ions that have 
been previously described in All and Be nson et al.| 
(2011). We also provide a csb measurement for SPT-CL 
.12106 —5844, whose o ther X-ray properties are discussed 



i 



Foley et al. (20111. Their csb values are provided in 



Table 

In the SPT sample, we find two high-redshift galaxy 
clusters with surface brightness concentrations in the 
strong cool-core regime. These are among the first strong 
cool-core galaxy clusters detected at redshifts beyond 
z = 0.5. One of the strong cool-core clusters is SPT- 
CL J2331— 5051, which lies at z — 0.58 and has a csb 
of 0.214 ± 0.016. The strongest cool-core cluster is SPT- 
CL J0000-5748 and has a c SB of 0.244 ± 0.023. At a 
redshift of 0.702, corresponding to 7.4 Gyr after the Big 
Bang, this cluster is also the highest redshift strong cool- 
core cluster found in this work. These two clusters rule 
out the hypothesis that there are no galaxy clusters at 
z > 0.5 classified by csb as having strong cool cores at 
the 5.4ct level. 

The X-ray properties of the gala xy cluster SPT-CL 
J2106-5844 were studied in detail in|Fole y et aT] (|2011|). 

10 ib h\ 



With a mass M- 



2i.hi 



(1.27 ± 0.21) 



T 
'70 



M: 



and a redshift z = 1.13, it is the most massive cluster 
known at z > 1. They measure a cluster temperature of 
Tx = H.Ol^g keV and a central temperature of Tx = 

6.5_ti'i keV within 0.17rsoo- The temperature decrement 
in the core suggests moderate cooling. We measure the 
surface brightness concentration csb = 0.026 ± 0.007, 
which lies below the moderate cool-core threshold. 

4.3. csb Distribution 

In this Section, we compare the csb distributions of 
the SPT and X-ray-selected hig h-redshift cluste r sam- 
ples. To remain consiste nt with Vikhlinin et al. (20071 
and Santos et al. (2010), we define a redshift break in 
our cluster samples at z = 0.5, with clusters at z > 0.5 
defined to be "high-z" and clusters at z < 0.5 defined to 
be "low-z". We do not include SPT-CL J2106-5844 in 
these analyses as it is not part of a uniformly SZ-selected 
sample with X-ray follow-up, as is approximately the case 
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Figure 3. Surface brightness concentration as a function of red- 
shift for objects with a measured c$b included in this study. Four 
galaxy clusters are members of two surveys and are shown with 
symbols from both surveys. Dashed line at cge = 0.075 corre- 
sponds to the boundary between noncool-core and moderate cool- 
core clusters. Dash-dotted line at csb = 0.155 corresponds to the 
boundary between moderate and strong cool-core clusters. 



for the clusters studied in All. 

The distribution of csb as a function of redshift for 
various X-ray-selected samples and the SPT sample is 
in Figure [31 At z > 0.5, there are 31 clusters in X- 
ray-selected samples and 9 in the SPT sample. In this 
regime, there are only two strong cool-core clusters, both 
in the SPT sample. This constrains the fraction of strong 
cool cores at high redshifts in the SPT sample to between 
7% and 56% at a 95% binomial confidence interval (see 
Camerori]|2011 for calculation details). 

The csb distributions of X-ray and SZ samples, broken 
up into low-z and high-z groups, are shown in Figure HJ 
We perform a K-S test between the SPT high-z sample 
and the X-ray high-z sample. The result yields a p-value 
of 0.42, indicating that although the SPT sample con- 
tains the only strong cool-core clusters, the two selection 
methods show no evidence for being drawn from different 
parent distributions. 

In order to study the evolution of the cool-core fraction 
we compare the csb distribution of the SPT high-z sam- 
ple with that of the X-ray-selected low-z sample (Fig- 
ure [4| 37 . A K-S test between the two datasets results in 
a p-value of 0.87, providing no evidence that these two 
distributions are drawn from different parent distribu- 
tions. Therefore, although there is a smaller fraction of 
strong cool-core clusters in the high-redshift SPT sam- 
ple, these results provide no evidence for evolution of the 
cool-core fraction between the two redshift regimes. 



4.4. Investigation of AGN Contamination 

Potentially, strong radio emission from AGN in cool- 
core clusters could affect their SZ detections. Previous 
investigations suggest that this should not be a signifi- 
cant concern for samples selected at 150 GHz (Lin et al. 
Sehgal et al.||2010 ). All investigate the role radio 



2009 



37 We do not study the evolution solely within the SPT sample 
because at z < 0.5, the sample has only four clusters. 
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Table 3 

SPT Galaxy Cluster Measurements 



Cluster 




R.A. 

(1) 


Dec 
(1) 


CSB 


z 

(2) 


BCG Offset 
(3) (kpc) 


SPT-CL JOOOO- 


-5748 


0.250 


-57.810 


0.244 


± 


0.023 


0.70 


7.4 


± 


3.5 


SPT-CL J0509- 


-5342 


77.338 


-53.703 


0.109 


± 


0.010 


0.46 


24.1 


± 


2.9 


SPT-CL J0516- 


-5430 


79.148 


-54.506 


0.026 


± 


0.005 


0.30 


107.3 


± 


2.2 


SPT-CL J0528- 


-5300 


82.020 


-52.997 


0.061 


± 


0.009 


0.77 


58.1 


± 


3.7 


SPT-CL J0533- 


-5005 


83.405 


-50.098 


0.013 


± 


0.007 


0.88 


414.2 


± 


3.8 


SPT-CL J0546- 


-5345 


86.655 


-53.759 


0.072 


± 


0.010 


1.07 


40.8 


± 


4.0 


SPT-CL J0551- 


-5709 


87.893 


-57.144 


0.034 


± 


0.008 


0.42 


82.0 


± 


2.7 


SPT-CL J2106- 


-5844 a 


316.518 


-58.742 


0.027 


± 


0.007 


1.13 


24.3 


± 


4.1 


SPT-CL J2331- 


-5051 


352.963 


-50.865 


0.214 


± 


0.016 


0.58 


4.5 


± 


3.2 


SPT-CL J2337- 


-5942 


354.353 


-59.705 


0.033 


± 


0.006 


0.78 


199.4 


± 


3.7 


SPT-CL J2341- 


-5119 


355.301 


-51.329 


0.092 


± 


0.009 


1.00 


3.0 


± 


4.0 


SPT-CL J2342- 


-5411 


355.692 


-54.185 


0.138 


± 


0.013 


1.08 


8.6 


± 


4.1 


SPT-CL J2355- 


-5056 


358.948 


-50.928 


0.113 


± 


0.010 


0.32 


6.7 


± 


2.3 


SPT-CL J2359- 


-5009 


359.931 


-50.170 


0.035 


± 


0.007 


0.78 


85.5 


± 


3.7 



Not e. — (1) Coordinates determined from X-ray centroid. (2) Rcdshifts fromlSongl 
|et al.| l [2012[l . (3) Pr ojected offse t between BCG position and X-ray centroid, with BUTT 
positions taken from pong et ah| | |2012| . The error in offset corresponds to the resolution 
limit of Chandra at the cluster s redshift. 

a SPT-CL J2106— 5844 is not included in the analysis of the csb distribution. 
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Figure 4. Distribution of the surface brightness concentration 
of the X-ray high-z sample (the z > 0.5 clusters in the CCCP 
high-z sample and the z > 0.6 clusters in the RDCS+WARPS 
samples; green forward-hatched areas), the X-ray low-z sample (the 
z < 0.5 clusters from the CCCP high-z sample and the CCCP low-z 
sample; red empty areas), and the SPT high-z clusters (the z > 0.5 
clusters from the SPT sample; blue back-hatched areas). Insert: 
Cumulative distribution of the surface brightness concentration for 
the same three samples. 



AGNs play specifically for the SPT samp le. Using the 
Sydn ey University Molonglo Sky Survey ( |Mauch et al. 
120031) at 843 MHz, All identify radio sources in the cen- 
ters of 7 of the 13 clusters from the SPT sample. They 
extrapolate the radio flux to the SPT bands and find the 
total flux of these sources to be negligible compared to 
the SZ signal. 

Measurement of csb can be affected by the presence of 
an X-ray-loud AGN within the core radius of the clus- 
ter. Such systems result in inflated csb values due to the 
AGN's contribution to the core X-ray counts. In this Sec- 
tion, we investigate AGN contamination in the moderate 
and strong cool-core clusters from the SPT sample. 

While there are many signatures that a galaxy con- 
tains an AGN, there is no definitive diagnostic that the 



AGN is emitting X-rays besides direct observation of X- 
rays. Therefore, with the high resolution data provided 
by Chandra, the best way of identifying such features is 
to visually inspect the images to investigate whether any 
point sources are within the inner 40 kpc of the cluster 
center. To validate the visual inspection, we also em- 
ploy a fully-automated source-detection method using 
the CIA0 tool wavdetect. We also specifically investi- 
gate for the presence of a point source coincident with 
the X-ray centroid by comparing the hardness ratios of 
the cluster center with an annulus around the center. 
The centrally decreasing temperature of a cool-core clus- 
ter means that an uncontaminated core will have softer 
emission in the center than in the annulus. An X-ray 
AGN, however, has hard X-ray emission and its presence 
in the cluster core will cause the core to have a higher 
hardness ratio than the annulus. 



We utilize BEHR ( Park et al.||2006| ), a fully Bayesian 
approach to calculate hardness ratios that treats pho- 
ton counts as Poisson statistics with appropriate error 
propagation. We use an X-ray hardness ratio defined 
as (H — S)/(H + S) where H corresponds to the X-ray 
counts in the hard band (2 — 8 keV) and S is the X-ray 
counts in the soft band (0.5 — 2 keV). The hardness ratio 
of the core is taken in a 2 arcsecond aperture centered 
on the X-ray centroid. This aperture size is used as it 
would capture approximately 90% of the X-rays from a 
point source located in its center. The annulus has an 
inner radius of 2 arcseconds and an outer radius of 4 
arcseconds. 

In Table |1J we give the results of our checks for AGN 
contamination. Visual inspection of the Chandra im- 
ages reveals no evidence of an X-ray point source within 
40 kpc of any of the cluster centers. The source-detection 
tool wavdetect only detected a source in the center of 
one cluster: SPT-CL J0000-5748. In this case, the ratio 
of the size of the detection with the size of Chandra's 
point spread function at the source location is 1.8, indi- 
cating the detected source is more extended than a point 
source (e.g., an AGN) and the emission comes from the 
cluster's cool-core ICM. In addition, the core hardness 



ratios range from — 0.73l ; 2 3 to 



-0.48±g;{£ and are all 
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Table 4 

AGN Signatures for Moderate and Strong Cool-Core Clusters 
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Cluster 




HRannulus Central Source 






Wav Vis 




(i) 


(1) (2) (3) 



SPT-CL J0000- 


5748 




-o.39±° : gs 


Yes 


No 


SPT-CL J0509- 


5342 


-0 51 + - 15 
u ' ol -0.15 


-0.08t ;l 2 2 


No 


No 


SPT-CL J2331- 


5051 


-0 58+° 11 


-0 V)+ 012 


No 


No 


SPT-CL J2341- 


5119 


-0 fi4+ 016 


-0 3^+ 012 


No 


No 


SPT-CL J2342- 


5411 


-0 fifi +0 13 
u - oo -0.19 


-0 >53+ 014 


No 


No 


SPT-CL J2355- 


5056 


_ n -O+0.13 


-0 52+ 018 


No 


No 



Note. — (1) Hardness ratio calculated in the cluster core and its 
surrounding annulus. A core with a lower hardness ratio than its annulus 
is consistent with no contamination from an X-ray AGN. (2) Indication 
of a source detection by wavdetect at the cluster center. (3) Indication 
of an X-ray point source within the inner 40 kpc of the cluster center 
based on visual inspection of the observations. 

lower than the hardness ratios in their surrounding an- 
nulus. These core hardness ratios are also softer t han the 
hardness rat ios of typical AGNs. For example, Hickox 
et al. (2009) report a hardness ratio of —0.37 ± 0.06 for 
a stack of 95 radio-selected AGNs. Therefore, we find 
no evidence for AGN contamination in our csb measure- 
ments for the SPT sample. 



5. DYNAMICAL STATE AND COOL-CORE STRENGTHS 

In this Section, we examine the relationship between 
a cluster's dynamical state and its cool-core strength. 
This can be used to investigate the effect of mergers on 
the cool-core properties of our sample. For this analysis, 
we characterize a cluster's dynamical state based on the 
offset of its BCG from the center of its X-ray emission 
( Katayama et al. 2003 ) . Having a complete sample is less 
important here because the BCG position is not involved 
in our selection techniques. Therefore, in this Section, we 
include SPT-CL J2106-5844 in the SPT sample. 

The BCG offsets are determined from the projected 
distance between the X-ray cen troid (as described for 



the csb cal culation in Section 3.1 ) and the BCG positions 
reported in Song et aL] ([2012 ). Uncertainties in the offsets 
are given as the resolution tor the Chandra observations 
at each cluster's redshift. 

Although the majority of local clusters exhibit rela- 
tively smal l BCG offsets ( |Lin fc Mohr1|2004[ ) , [Sanderson 



et al. ( 2009) ) show that the magnitude of this offset is an 
ticorrelated with the cool-core strength of the cluster. As 
shown in Figure [5[ the offsets for our SZ-selected sample 
are in agreement with these results. We find a strong 
anticorrelation between csb and the BCG offset as well 
as an absence of noncool cores with low offsets. The two 
parameters exhibit a Spearman's rank correlation coeffi- 
cient of p = —0.83 with a p-value of 2.2 x 10~ 4 . These 
results reinforce the model that as clusters relax, BCGs 
settle in the center of the potential well and cool cores 
become established. Later, dynamical disturbances are 
capable of removing a cool core. Within this model, the 
clean anticorrelation reinforces csb as a parameter that 
accurately measures cluster cores at high redshifts. The 
lack of "relaxed," noncool-core clusters with small BCG 
offsets suggests that neither prehe ating (as proposed by 
Kaiser|199l[ McCarthy et al.|2008|) nor strong AGN feed- 
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Figure 5. Offset between the X-ray centroid and the BCG as 
a function of cge- Errors in the BCG offset correspond to the 
resolution limit of Chandra at the cluster's redshift. Data from the 
SPT sample and SPT-CL J2106-5844. 



back (as reviewed in Fabian 2012; McNamara & Nulsen 

II r 



2012 1 are the most dominant mechanisms of cool-core 
disruption in our sample. Rather, this anticorrelation 
points to recent mergers playing a major role by mixing 
or shock-heating the ICM. 

6. CONCLUSIONS 

We study the cool-core characteristics of an SZ-selected 
cluster sample detected by the SPT survey. The sam- 
ple provides an alternative to the existing X-ray-selected 
samples for studying the cool-core fraction at high red- 
shifts, which may be impacted by X-ray selection effects. 

We evaluate the two parameters that are suitable for 
measuring cool cores at high redshifts, the concentration 
of surface brightness (csb) and the cuspiness (a), and 
examine their performance as a predictor of the central 
cooling time i CO oi- A Spearman's rank test shows correla- 
tion coefficients between a and t coo \ and between csb and 
icooi to be p a = —0.88 and p CSB = —0.84, respectively. 
However, csb exhibits much smaller fractional measure- 
ment errors. 

Using csb, we find evidence of two strong cool-core 
clusters in the SPT sample at z > 0.5, among the first 
of their kind. We rule out the hypothesis that there 
are no such galaxy clusters at the 5.4a level. With a 
sample of 9 high-redshift clusters, we show that the high- 
z strong cool-core fraction is greater than 7% with 95% 
confidence. We compare the distributions of the SPT 
sample and previous X-ray-selected samples with a K- 
S test. The result yields a p-value of 0.42, providing no 
statistically significant evidence that the two samples are 
drawn from different distributions. 

We also evaluate the relationship between the strength 
of a cluster's cool core and its dynamical state. We find a 
strong anticorrelation between csb and the offset between 
the BCG and the X-ray centroid, which is related to the 
dynamical state of the cluster. While preheating or AGN 
feedback may be responsible for the lack of cool cores in 
some clusters, this result suggests that in our sample, the 
formation of cool cores is inhibited by merger events in 
clusters, which cause turbulent mixing or shock-heating 
of the ICM. 
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The results of this study are based on a complete 
sample of SZ-selected clusters with Chandra observa- 
tions sampled from the first 178 deg 2 of the SPT survey. 
The statistics of our results will soon be dramatically 
improved through a Chandra X-ray Visionary Project 
underway to observe the 80 most significant clusters at 
z > 0.4 from the first 2000 deg 2 of the SPT survey. 
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